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REVISION OF GEODETIC PARAMETERS 
E. M. Gaposchkin and M. R. Williamson 

ABSTRACT 

Laser data from nine satellites and 12 stations are combined with surface- 
gravity data to ot In spnerical harmonics representing the geopotentlal com- 
plete through degree ind order 18. This laser-data-only solution provides a 
reasonable Improvement to the gravity field. 


INTRODUCTION 


Smithsonian Astrophyslcal Observatory has published gravity-field solutions 
utilizing both satellite-tracking and terrestrial gravity data (see, e.g., 
Gaposchkin, 1974), which were based primarily on precision-reduced camera data 
and on then-available surface-gravity data. Now, however, we have better data, 
and new types of data will soon be available. The work reported here Is the 
beginning of a general revision and extension of our knowledge ot the geopoten- 
tial. In this first iteration, we will experiment with new data (laser ranges), 
verify new methods of data reduction, and prepare for new types of data (altim- 
eter and satellite-to-satellite tracking data). We seek here to Improve our 
knowledge of the gravity field so that better satellite orbits, consistent 
with surface-gravity data, can be calculated. This is done by using laser data 
only, whose accuracy (approximately 1 m) is far greater than the ephemerls 
accuracy (approximately 10 m). The situation has changed since 1971, when the 
bulk of the data used were 4-arcsec camera data. Furthermore, substantial re- 
visions have been made in the treatment of orbit perturbations, and all these 
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advances indicate a new attempt at improved geopotential coefficients. 

OBJECTIVES 

This computation has several objectives, the primary one being to use 
laser data only in a determination of the earth's gravity field, with the aim 
of computing satellite orbits to an accuracy comparable with that of the laser 
data. We can obtain a realistic gravity field consistent with surface-gravity 
data. A secondary objective is to study the consequences of using data that 
provide no ties to an inertial reference frame, as was the case with camera 
data. For our third objective, we will investigate the effects of one satellite 
on another in the solution; that is, we can optimize a solution for one 
satellite by using only data from that satellite. The question 
then becomes: How much does adding data from a second satellite degrade the 
orbit computed for the first? Of course, improved orbits and a more accurate 
geoid are necessary for analyzing satellite altimetry data for geodetic and 
oceanographic purposes. 


REFINEMENTS AND TECHNIQUES 

Improvements in perturbation calculations have been numerous. The inclination 
function for tesseral harmonics, as formulated by Kaula, computationally loses 
accuracy for high-degree coefficients. It has been replaced by a mathematical- 
ly equivalent formula derived from group theory (Gaposc : >kin, 1973). The inter- 
action terms between and resonant harmonics have also been improved. Lunar 
and solar perturbations and body tides and ocean tides have been computed to the 
necessary accuracy (Kozai, 1973). Perturbations arising from the noninertlalness 
of the coordinate system have been corrected and improved (Kinoshita, 1975a, b). 
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and those due to direct sol.'r radiation pressure (Aksnes, 1975), albedo 
pressure (Lautman, 1975a, b), and Infrared radiation (Lautian, 1975c) have all 
been Included and tested (Gaposchkin et al . , 1975). 

The compilation of surface-gravity data used in Gaposchkin (1974) Pas since 
been augmented (Williamson and Gaposchkin, 1975). The surface-gravity data are 
summarized in Table 1, and the distribution of these data is shown in Figure 1. 
Table 2 compares the 1° x 1° gravity anomalies from the Defense Mapping Agency 
Aerospace Center (DMAAC, 1973) with other available compilations. 

Three coordinated programs have provided laser tracking data in sufficient 
density to be used for a gravity-field determination: the International 

Satellite Geodesy Experiment (ISAGEX) held in 1971, the Earth Physics Satellite 
Observation Campaign (EPSOC) in 1972-i973, and the Geos 3 campaign in 1975. 

There has been a steady improvement in the volume, reliability, and accuracy of 
the data. One of the objectives of the ISAGEX program was to obtain data for 
determining the gravity field, and a number of orbital arcs are suitable for this 
purpose. The EPSOC program was directed to the study of long-period effects and 
polar motion. Some arcs from EPSOC are also used for determining the tesseral 
harmonics. Finally, during the current Geos 3 program routine data are being 
obtained and Included in the analysis. Table 3 lists the participating stations 
for each campaign that are used in the analysis reported here. The distribution 
of the stations is shown in Figure 2. 

Currently, nine satellites in orbit are equipped with cube-corner reflectors 
suitable for laser ranging; these satellites are useful for a gravity-field de- 
termination, especially as their distribution in inclination and height is 
reasonably good. Table 4 list 1 : their orbital characteristics and the number of 
arcs used, and the distribution of satellites is given in Figure 3. 


Satellite orbits are usually computed for between 8 and 16 days, the In- 
terval depending primarily on the availability of data. Also, each arc covers 
at least one oscillation of the resonant period. Table 5 gives the resonant 
periods for the satellites, while Table 6 lists the constants used to calculate 
the satellite orbits. 


SOLUTIONS 


The normal equations for surface-gravity data have been computed complete 
from degree 2 through 18. The combination solution Included a number of harmon- 
ics of higher degree that are resonant with one or more of the satellite orbits. 

To this set of surface-gravity normal equations was added a set of normal equa- 
tions, satellite by satellite. The system was solved after adding one satellite 
(Geos 1), and the result has been compared with all the satellites and with 
surface gravity. For the satellite comparison, one arc from each satellite 
was selected: generally, the arc with the most data. The orbit was recalcula- 

ted with the revised gravity field, and the orbital fit in terms of o Q was used 
as the criterion. In all combinations, the satellite data were used at their 
a priori weight. The surface gravity was given se eral weights, and all solutions 
were tested in order to determine the optimum weight for the combination. 

The weight finally adopted is 
27 

< A J nA mgal ’ 

where n is the number of 1° x 1° squares in each 5° x 5° mean, A is the area 
of the anomaly, and <A> is the average area. This weight is twice that used in 

the 1973 Smithsonian Standard Earth (III) (SE III) (Gaposchkin, 1973). The 

ISAGEX laser data were given a 5-m weight, and all other laser data, a 2-m weight. 

These combination solutions are si'.nmari zed in Tab’ 7. The orbit for a satellite 

not used in the solution is really ver> poor, generally because of relatively small 

changes in a few coefficients resonant with the gravity field. We note that 
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Peole, which has no resonance, degrades marginally. The surface gravity 
<(Ag - Ag ) 2 > Is only relative. A large part of those residuals Is Information 

In the higher harmonics. The best fit Is obtained with only one satellite; 

however, when satellite data are added, the degradation Is not large considering 


the overall accuracy of surface-gravity data. 

Generally, adding satellite data also degrades the satellite orbit fit, but 

not very much. This overall Improvement is considered quite satisfactory for 

one iteration of a very complex nonlinear process. All satellite orbits Improved 
2 

by at least 1 m In the orbital fit. In percentage terms. Geos 1, BE-C, Geos 2, 


Starlefte, Peole, and Geos 3 each improved in orbital fit by 17.9, 22.2, 43.9, 


69.2, 26.2, and 13.6%, respectively, for an average 32% improvement! 

The final adopted 52-arc solution (SE IV. 1) can be compared with the surface-gravity 
data In more detail. Assuming they are statistically independent, the following 
quantities defined by Kaula (1966) can be computed and use. to compare a geopo- 


tential model (y^) with observed values of surface gravity (g^): 

(g 2 ) The mean value of gj 2 , where g t is the mean free-air 

gravity anomaly based on surface gravity, indicating the 
amount of information contained in the surface-gravity 
anomalies. 
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The mean value of g , where g is the mean free-air 

6 S 

gravity anomaly computed from the gcapotential model, 
indicating the amount of information in the computed 
gravity anomalies. 


<g t S s > 


<(g t - g s ) 2 > 

E(e 2 ) 

E (6g 2 ) 


An estimate of g h - i. e. , the true value of the < ontributicn 
to the gravity anomaly of the geopotential model and the 
amount of information common to both g^ and g g . 

The mean-square difference of g t and g R . 

The mean-square error in the geopotential model. 

The mean-square error of the observed gravity. 

The mean square of the error of omission - that is, the 
difference between true gravity and g^; this term is then 
the model error. 
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If the geopotential model were perfect, then <g?> ■ <gjj>, wh 1ch In turn 
would equal < 9 t 9 s > If 9 t were free from error and known everywhere. Then, 
c j would be zero even though g $ would not contain all the Information necessary 
to describe the total field. The information not contained In the model field — 
l.e., the error of omission, 6g - then consists of the hlgl er order coefficients. 
The quantity <(g. - g ) 2 > is a measure of the agreement between the two estimates 
g t and g $ and Is equal to 

«g t - g s ) 2 >= E(^) ♦ E(^) 4 E(6g 2 ) . 

Another estimate of g^ can be obtained from the gravimetric estimates of 
2 

degree variance (Kaula, 1965 ); 

E( *h>- D -Zrrr«? . 

/ 


where n c is the number of coefficients of degree t Included in g^, and 


2 o 
V ('-!)“ 


Eft. ♦■?».) 


m 


We also have 

E(e s> 51 <R 2 s ) ~ <g g g t ) 

and 

E(cf) = <«f>/<n) . 


These values are given in Table 8 for SE III and for this solution. 
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The Information In the surface-gravity data solution <g|> has Increased 
In this new data set. This Is reasonable; since the unobserved areas had an 
expected value of zero, the fewer observed areas there are, the lower the 
variance Is. However, the Information In the satellite solution <g£> has 
decreased, a fact that Is confirmed by a decrease In D. Therefore, the Infor- 
mation In SE III was too high. The residual <(g t - g s ) 2 > has remained roughly 
the same, while the Information In the higher harmonics Is estimated to be 
’argei. Ire estimate of E(e^) cannot be good, as these sets of data, g s and 
g t , are not independent. 

The spherical-harmonic coefficients are listed In Tables 9 and 10. Flqure 4 
is a plot of the mean potential coefficient as a function of degree, and 
Figures 5 and 6 show the geoid height and gravity anomalies for this solution. 

FUTURE WORK 

The obvious next step is to complete another Iteration, taking the solution 
to perhaps degree and order 24. To this can be added the normal equations for 
zonal harmonics and sets of resonant harmonics. When the orbital accuracy ap- 
proaches a few meters, then we must reduce the error In the station coordinates by 
solving for them. Finally, of course, the aim is to add altimetry data to this 
system of normal equations. 

In summary, we find that the improved accuracy that has been apparent from 
laser data is becoming realized. 
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Table 1. Surface-gravity data available. 


Source 

1° 

No. of 
x 1° means 

ttO. flf. S5Q E 

n 2 1 

km 

n ■ 25 

Gaposchkln (1974) 

19328 

1183 

145 

Williamson and Gaposcbkln 

(1975) 

31636 

1452 

485 

Maximum number 

64800 

1690 

1690 
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Table 2. Comparison of 1° x 1° mean gravity anomalies with DMAAC (1973). 


Source 

No. of points 
compared 

Mean difference 
(mgal) 

rms 

(mgal) 

Australia (Mather, 1970) 

1364 

1.64 

24.16 

North America and North Atlantic 

(Talwanl et al.. 1972) 

3613 

-0.18 

15.29 

Indian Ocean 

(Kahle and Talwanl, 1973) 

2226 

-1.66 

23.09 

Worldwide (ACIC, 1971) 

19164 

-0.23 

16.99 
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Table 3. Lasers stations 

used in 

this analysis. 



. Station 


I SAGE X 

EPSOC 

Geos 3 

Number 

Location 

Agency" 




7902 

Ollfantsfonteln, S. Africa 

SA0 

X 

X 

X 

7907 

Arequlpa, Peru 

SA0 

X 

X 

X 

7921 

Mt. Hopkins, Arizona 

SAO 

X 

X 

X 

7928 

Natal , Brazil 

SAO 

X 

X 

X 

7930 

Athens, Greece 

SAO 

X 


X 

7050 

GSFC 

NASA 

X 


X 

7060 

Guam 

NASA 

X 



7061 

San Diego, Calif. 

NASA 


X 


7080 

Quincy, Calif. 

NASA 


X 


7068 

Grand Turk Island 

NASA 



X 

7804 

San Fernando. Spain 

CNES 

X 



7809 

Haute Prov • . • France 

CNES 

X 
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Table 4. Summary of dynamical data. 



Satellite 

Incl Inatlon 

Eccentricity Perigee 
height 
(km) 

(km) 

Number 

of 

arcs 

Designation Name 

7010901 

Peole 

15° 

0.017 

635 

7070 

5 

6701401 

DID 

39 

0.053 

569 

7337 

3 

6701101 

D1C 

40 

0.052 

579 

7336 

2 

6503201 

BE-C 

41 

0.026 

941 

7311 

9 

7501001 

Starlette 

50 

0.0207 

805 

7335 

5 

6508901 

Geos 1 

59 

0.073 

1121 

8074 

14 

750C701 

Geos 2 

-65 

0.0005 

840 

7222 

4 

6800201 

Geos 3 

-75 

0.031 

1101 

7709 

8 

6406401 

BE-B 

80 

0.012 

912 

7362 

2 






Total 

52 
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Table 5. Resonant periods 


Satellite 

Resonant with 
order m 

Period 

(days) 

7010901 

none 


6701401 

13 

9.4 

6701101 

14 

2.6 

6503201 

13 

5.6 

7501001 

14 

3.2 

6508901 

12 

7.2 

7502701 

14 

3.9 

6800201 

13 

6.3 

6406401 

14 

2.9 


l 
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Table 6. Constants used In orbit computation. 

GM - 3.986013 x 10 20 cm 3 sec -2 
c * 2.997925 x 10 10 cm sec' 1 * speed of light 
* 0.25 * Love's number 
c£ * 10° ■ phase lag of tide 

a e « 6.378140 Mn 
a * 0.32 ■ earth's albedo 


Satellite 

G 

«o 3 
1 

7010901 

0.20 

6701401 

0.30 

6701101 

0.30 

6503201 

0.13 

7501001 

0.01 

6:089Pi 

0.10 

7502/01 

0.04 

6800201 

0.06 

6406401 

0.10 
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Table 7. Comparison of solutions. 





I 




oo 

vO 

** 

VO 

o> 

00 

OJ 


** 

ro 

# 

rH 

VO 

vO 


C 


ro 

ro 

ro 

CM 

CM 

4~* 

07 


CM 

CM 

CM 

CM 

CM 

CM 



rH 



|H 

VO 

CM 

»T3 


vO 



rH 

Ov 

eH 

• 

UO 

o 

O 

VO 

»H 

CM 

770 

oo 

rH 


VO 






CM 

VO 

o 

cO 

o 

CM 





rH 

CM 

**■ 

CM 

ro 

CM 

CM 



C 

rH 

rH 

rH 

rH 

rH 

•"H 





rH 

rH 

rH 

rH 

rH 


on 

r\ 










<5! 




00 

o 


ro 

CM 

OO 



o 

vO 

vO 


CM 

ro 





o 

t 

• 


• 

• 

. 





CO 

CM 

ro 

ro 

oo 

ro 






ro 

oo 










rH 








o> 

8 


CTv 

VO 

to 





o% 

o 

oo 

o> 

CT» 





VO 

vO 


V£> 

VO 

VO 





rH 

•"“* 

rH 

rH 

rH 

rH 

1 

tU 










CO 




rH 

r^. 

CM 

O 

o 

CM 




o 

oo 

00 

CM 

CM 

VO 

vO 




o 

• 

• 

• 

• 

• 

# 





CO 

00 










oo 










ro 


CM 

rH 

CM 





CM 

ro 

§ 

rH 

rH 

CM 

• — « 



c 

o 

O 

o 

o 

o 





ro 

ro 

ro 

ro 

ro 

ro 

<A 










>— p 




ro 

CO 

CM 

00 

00 

CM 



o 

CO 

rH 

rH 

ro . 

uo 

r* 




o 

• 

• 

• 

• 

• 

• 


00 

rH 

v| 




ro 

oo 

ro 

oo 

oo 

>> 









4-» 

<W 


o 







•r- 



CM 







> 

03 

m 

A 


A| 

CO 

ro 

8 

8 

CT* 

66 

§ 

W <N 

CM 

rH 





rH 


r" 

c 








VA 

*3 








07 

o> 

o> 








U 

< 

E 








*3 




ro 

ro 

O 

CM 

o> 

CM 

»♦- 

1 


A | 

0> 

CM 

ro 

ro 

ro 


U 




rH 

r— H 

rH 

rH 

rH 

rH 

=3 

4-> 


c 







CO 

cn 










< 




















V 









o 

v/> 



ro 

VO 

CO 

CO 

'fcioocvimMinui 



u 



O 






• 

u 



CM 






o 

to 









7ZL 








07 


c 








4~> 


o 








4-> 


' r— 




•— • 

rH 


CM 

rH CM 07 

ro 

4-> 




rH 




07 »— 


3 




i i 

VA 

u 

<A 

VA U (A CO r— W 

VA 






O 

T3 | 

TJ O 

O 1 O 1 O <J O <T3 

X3 O 

O 




UJ 

5 

C UJ 

C 01 

07 UJ 07 UJ rH rH 07 4-> 
UCQOCQQOOl/) 

C 07 

CO 




CO 

to co 

n o 

*3 O 


VO 



Table 8. Comparison of surface gravity -ith solution (-9-1*); all solutions are to the 18th degree and order 



Table 9. SE IV. 1 fully normalized tesseral-harmonic coefficients for the geopotential (in units of 10 
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FIGURE CAPTIONS 


Figure 1. Distribution of 1° x 1° mean surface-gravl ty data. 

Figure 2. Locations of the observing stations Included In SE IV. 1. 

Figure 3. Distribution of perigee heights and Inclinations of the satellites 
used In SE IV. 1 . 

Figure 4. Mean po.ential coefficient by degree. 

Figure 5. SE IV. 1 geoid height In meters calculated with respect to the best- 
fitting ellipsoid, f * 1/298.256. 

Figure 6. SE IV. 1 gravity anomalies in mllligals calculated with respect to 
the best-fitting ellipsoid, f * 1/298.256. 
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Figure 1. Distriiwtion of 1° V l e mean surface-gravity data. 
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Figure 5. SE IV. 1 geoid height in meters calculated with respect to the 
best-fitting ellipsoid, f = 1/298.256 




Figure 6. SE IV. 1 gravity anomalies in mi 1 1 i gals calculated with respect 
to tno best- fitting ellipsoid, f = 1/298.256 



